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!>■ ■ Abstract 

On 

The interaction due to the exchange of massless neutrinos between neutrons 
Q_|' is a long-range force. Border effects on this multibody exchange inside a 

dense core are studied and computed analytically in 1 + 1 dimensions. We 
demonstrate in this work that a proper treatment of the star's border effect 
automatically incorporates the condensate contribution as a consequence of 
the appropriate boundary conditions for the neutrino Feynman propagator 
inside the star. The total multibody exchange contribution is infrared-safe 
and vanishes exactly in 1 + 1 dimensions. The general conclusion of this work 
is that the border effect does not modify the result that neutrino exchange is 
infrared-safe. This toy model prepares the ground and gives the tools for the 
study of the realistic 3 + 1 star. 
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The possible connection between the multibody neutrino exchange and the stability of 
compact stellar objects has recently motivated several relevant works 0]- [||]. The interest 
raised by Fischbach's original idea JTJ is based on the presumably catastrophic consequences 
for the self-energy of compact objects, such as neutron stars, originated by the long-range 
neutron correlation due to massless neutrino exchange. This catastrophic effect is then 
invoked to justify the introduction of a lower bound for the neutrino mass. In a previous work 
0, we have shown that the total contribution of the many-body massless neutrino exchanges 
may be directly and rather easily computed, using an effective Lagrangian, and results in an 
infrared well-behaved star self-energy The catastrophic result of the resummation method 
is simply due, in our opinion, to the fact that it is done outside the radius of convergence of 
the perturbative series. 

Smirnov and Vissani ||, following Fischbach's approach of summing up many-body 
exchange, order by order, showed that the 2-body potential is damped by the blocking effects 
of the neutrino sea || and conjectured that a similar effect for a many-body potential would 
reduce Fischbach's catastrophic effect. 

The effect of such a condensate has also been incorporated in our non-perturbative 
method by using a neutrino Feynman propagator inside a dense stellar medium with a 
condensate term M. However, this condensate does not bring any major modification to 
our conclusion that the total result of the multi-body massless neutrino exchange is infrared 
well-behaved. 

The objection to our previous work |JJ is that we had worked in the approximation where 
the presence of the neutron star border was negligible. In fact, we stressed in || our belief 
that the neutrino condensate had to be understood as a manifestation of the star's border, 
and a preliminary proof of that statement can be found in ref. [B]. 

In this paper, we shall demonstrate that a proper treatment of the star's border effect 
automatically incorporates the condensate contribution. And that this is a consequence of 
the appropriate boundary conditions for the Feynman propagator of the massless neutrino 
in the interacting medium. Finally, we will show that the border effect does not modify our 
main result that neutrino exchange is infrared- safe. 

Our tool in this proof is the computation of the Feynman propagator for the massless 
neutrino in the neutron star medium. In ||, ignoring the border effect and hence using the 
translational symmetry, we have found the following effective propagator 

(1) 



where try accounts for Z°-exchange diagrams between the neutrino and the neutrons of the 
medium in the static limit. 

The latter is an infrared regular propagator, which can be written as follows: 



<j{ — by q\ + ise 

where (g>) M = (go — b,q). The term ise is the infrared regulator, where s is a sign that 
stands for the appropriate "time convention" . The determination of s needs some care: the 
propagator has been rewritten in terms of the four-momentum g>, its "time" component 
being (q — b) instead of q . As we will show, the appropriate boundary conditions, i.e. the 
"time convention" , should be imposed to keep the usual distribution in the go complex plane 
of the poles for the Feynman propagator: since neutrino (antineutrino) states correspond 
to the positive energy (the hole of negative energy) solutions, then the positive (negative) 
energy propagates forward (backward) in time. This rule implies 

s = sign(g )sign(g - b) . (3) 

It is crucial to insist that by energy we here meanstrictly go an d not the combination (go — b), 
which expresses the distance of the energy level to the bottom of the potential V = b. 
As a consequence, b being negative^], we get these dispersion relations for neutrinos and 
antineutrinos: 

for \q\ > \b\ 

V|gl or, (4) 

for |g| < |6| 

Finally, the poles are in the first and/or third quadrants in g complex plane, as shown 
in fig. [3]. Had we taken s — 1, for a range of values of \q\ < \b\ poles would have appeared in 
the other quadrants (see fig. p. Such a pole should be taken into account when rotating to 
Euclidean time. This has not been done correctly in 0, inducing a minor error which will 
be discussed below. 

Let us now return to the central issue of this note, which is the estimation of the border 
effect. As an attempt to have a simple analytic result, we start, in this work, our study 
with the 1 + 1 dimensions toy model for the following reasons: i) we will be able to work 
analytically to the end, ii) the subtle problem of the ie convention is easier to master, 
iii) thanks to the relative simplicity of the problem, we will be able to use different comple- 
mentary techniques and to get a deeper understanding of the physics, which will prove most 
helpful for a more realistic problem, i.e. 3 + 1 dimensions fTO . 




1 In ||, eq. (2), uj = \q\ ±6 @, is misleading. In fact, in an infinite star, the choice between applying 
Feynman's prescription to go or (go — b) might appear as open, but when boundary conditions are 
properly taken into account, the choice of go becomes mandatory. 
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FIG. 1. Schematic distribution of the poles of the propagator. The black dots and the grey 
squares represent the poles placed by following the time convention ise and +ie, respectively. 

For simplicity, we consider a sharp border located at z = and use an effective neutrino 
Lagrangian that summarizes the interaction with the neutrons H: 



£ e jf = iis(x)ftv(x) — b8(z)u(x)'y u(x) , 

where b ~ — ^IGp^f ~ —0.2 x 10~ 7 GeV, n n being the neutron density of the star 
The usual definition of the propagator is the following: 

S F (x,y) = (0\Th(xMy))\0) . 



(5) 



(6) 



It is worth stressing at this point that the time-ordered product will give us the Feynman 
boundary conditions. 

Following Gavela et al. f|, \l/(x), the general Dirac solution of the problem, can be 
written as: 



*(x) = £ b h (±k)v h ±J ;(x)e- iEt + di(±k)u^ h (x)e 



iEt 



(7) 



±k,i, 



where v are the eigenstates of the Dirac Hamiltonian derived from the Lagrangian (H). In the 
second quantization, the coefficients b and S become annihilation and creation operators, 
respectively. 

The functions v can be obtained from the following complete set of solutions: 



u^(z 7 t) = e- iE ^u h (±k) {9(-z) + e^6(z)} , 

u K h ( z , t) = e iE ^ z \(±k) {e(-z) + e^ bz 6(z)} , (8) 

which are obtained by solving the equation of motion for the 1 + 1 Lagrangian of eq. (H), 

{j°(id t -W(z))+i'y l d 1 }u(x,t) = Q . (9) 

In eq. (|8|) , the index h expresses the chirality of the solution and A denotes the antipar- 
ticles. Note that the chirality is the same on both sides of the border, as expected from the 
chirality-conserving Lagrangian of eq. @ . 

In the region z > (inside the star), neutrinos have k z = ±(E — b) and antineutrinos 
have k z = ±(E + b), while in the region z < (outside the star), they have both k z = ±E. 
The subscripts ±fc depend on the choice of the sign for k z . In Eq. (|8]) u and v are the Dirac 
spinors for particles and antiparticles, respectively, 

where fc> = (E — b, ±(E — b)) for neutrinos, fc> = (E + b, ±(E + b)) for antineutrinos and 
k% = (E, ±E) for both. 

Equations (|7|) and @ imply a definite choice of the zero energy level that corresponds 
to the standard choice of the free neutrinos far outside the star, z —>■ — oo, clearly the only 
admissible choice. Physically this choice is related to the use of a stationary Lagrangian, (0), 
which means that we assume that the system has relaxed to equilibrium, implying plane- 
wave solutions that extend outside and inside the star. These plane waves "know" the zero 
energy level from their outer domain. This choice, combined with the time-ordered product 
of eq. (H), implies the +ise time regulator in momentum space as announced before. 

In momentum space, the propagator can be written as 



S F (qf, g*)7° = i Jo dh \v h ±% {ql) {vK{q\ 



i 



j 

qo+ko—ie 



qo~h)+ie . . 



where q^ % = (go, Qz )■ Notice that the character of creation or annihilation for the operators 
in eq. ((?]) is fixed by taking the matter-free vacuum as reference, i.e. it is a consequence of 
our choice of the zero energy level. V are the Fourier transform of the eigenstates v{q z ) = 
f+^ dz e-">-^(z): 



We took E = ko in order to make the notation uniform. Applying this last result in eq. (|TlD , 
with the appropriate analytic continuation of the integration variable ko, and by following 
the adequate integration contour (see appendix |5|), we obtain for the 1 + 1 propagator: 



b i 



S F ( q ', q ')=2*S( qi - qh ^ , ,-j-^ 



where we define: 



tr^t (1 + sign ^°)^) + W^f< (1 ~ sign ^°)^)} 



r = i. = 2 i • with (s<) M = ( f = (so, g) 

?< <?< + *e 
1 <k 



( ^-^T-with (*>)" = (*-&,*) , (14) 

with s = sign ( g ) sign (g — b) and a z = 7°7 1 . The location of the poles in the complex plane 
have already been depicted in fig. |1|. 

Anticipating over the vacuum energy calculation, it is interesting to notice that the 
second term in the r.h.s of eq. ( |13|) does not contribute when the border is sent to infinity; 
one is then left with only the first term, which is exactly the one of eq. (0) for the infinite 
star with the same good regularization. To convince ourselves, eq. ( |13D can be rewritten as 
follows: 






(15) 



where V stands for the principal value. For an observer located at z — > oo, which means that 
he does not see the star border, restoring the translational symmetry, i.e. q{ — > q l z (when 
z — > oo), is a good approximation. Consequently, the second term in the r.h.s of eq. flT5|) 
will not contribute when z — ► oo 0. Furthermore, it can be seen from eq. (0) below that 
after the integration over the momenta, as suggested by Schwinger's method to obtain the 
vacuum energy, the second term in the r.h.s of eq. (|13D does not contribute to the energy 
density w(z — > +oo). 

The expression for l/(^)*, given by eq. (fll|), can be appropriately rewritten as 



2 In fact, we can see in eqs. (ID]) and (O) below that the break-up of the translational symmetry is 
expressed in e l ^ qz ~ qz ' z , which gives the dependence on z. This oscillating term, in the limit z — ► oo, 
will destroy any contribution to the integral coming from the part of the principal value of the 
propagator. 
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2iri4 > 6(<£)9(qo)0(b-q o )\ . (16) 

This last equation is a concrete demonstration that we have, in 1 + 1 dimensions, 
generated the condensate contribution in a natural way to all orders. In refs. || and [Bf, 



the same condensate term in the r.h.s of eq. (16) was introduced by hand. We should stress 



that in ref. we initially computed the weak energy density by using the +ie convention 
(first term of eq. (|16|)), but we did not take into account the pole go = \o\ + b, in the 
second quadrant for \q\ < \b\ (see fig. [I]) when rotating to Euclidean time. This explains our 
non-zero result for the energy density when we have fixed the good boundary conditions by 
introducing the second term of eq. fllTf ) for the propagator. 

The existence of a certain neutrino condensate due to the interaction of neutrinos and 
the stellar matter background was initially proposed by Loeb M; neutrinos are trapped 
inside the star and antineutrinos are repelled. The relevant physics concerning the neutrino 
propagator in the scenario we described above would appear by introducing the appropriate 
boundary conditions. Equations ( |il| ) and flT6"D give a confirmation of the idea proposed in 
ref. @ : the condensate is a consequence of the existence of a border. 

This condensate is physically understandable. As we have tuned the level of the Dirac 
sea outside the star (to the left) and as our states extend over all space, far inside the star 
(to the right), the level corresponds to filling a Fermi sea above the bottom of the potential 
b < q < 0. This obviously induces a Pauli blocking effect^. Equation ([|) anticipates this 
result: \b\ is a lower bound for the momentum of the go-P°shive states. 

With these tools, let us compute the total neutrino-exchange contribution to the energy 
of the star, which is nothing else than the difference between the vacuum energy of our 
effective theory eq. (^) and the vacuum energy of free neutrinos (without the star). Following 
Schwinger's method [0], the vacuum energy in the presence of an external field is given by 
tracing the Hamiltonian multiplied by the propagator, as done in refs. fl|] and ||. Now, we 
compute the time-independent vacuum energy density 

w(z) = -idMl° (S F (z, t; z>, t) - Sf{ Zl t; z>, t)) P L } Z ^ Z 

■ qol o e i(d- q i)z (s F ( q f } qi) _ 5-(o) (g/j qi) 2j[6 ( q f _ q ty Pl 



dq{ dq\ dq , 



(2tt)3 



(17) 



where S F is the free neutrino propagator and Pl = —^. We should remark that the 
e i{qz-ql)z j n th e second line of this equation is generated by the break-up of the translational 
symmetry due to the border. 



3 Smirnov and Vissani [fj] have proved that the condensate term, obtained above in a natural way, 
generates the damping of the 2-body neutrino-exchange potential. 



Now, by applying the expression of eq. (|1 
obtain: 



for the propagator in this last equation, we 




^\Pr 



1 + sign(g )a 2 



iiy 
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1 -sign(g )a 2 ) \ Pl 



(18) 



It should be noted that the first term of the r.h.s can be appropriately rewritten by using 
the Schwinger-Dyson relation, which has been shown in ref. |2| for the effective propagator 
(H) in the case of the infinite star []. However, the calculation is rather cumbersome, even if 
it presents no special difficulty and leads to the amazing result: 



w(z) = 



(19) 



As surprising as this result may seem, it can be understood physically in a rather simple 
way. One must remember that w(z) is nothing else than: 



t 



E n m{ (z)V, 



z - 



E^lHz)^_(z) 



(20) 



where n_ labels the negative energy states and (0) refers to the matter-free vacuum. 

These integrals are ultraviolet-divergent and we will allow the interchange of the sum 
and the integral as a regularization method. This was implicitly done in eq. fll8l) . Looking 
now at the solutions of eq. (Q) of the Hamiltonian with and without (b = 0) the star, we 
first that remark there is a one-to-one correspondance of the states in the two situations and 
with the same energy; furthermore the probability density of the states are all equal to 1 for 
all z. More explicitly, the border does not disturb the wave functions, except for a phase. 
Consequently, after regularization (interchanging the sum and the difference) each term in 
the sum of eq. (^) vanishes. 

The (1 + l)-dimensional problem presents the particularity that, for a given wave- 
plane solution, the sign of the momentum determines a positive correspondence between 
z = + (— )oo and t = + (— )oo, or vice versa. Massless fermions are not reflected by a 
one-dimensional potential. However, the potential introduced through the Lagrangian (|5]) 
does not allow us to identify, for a given process, asymptotically free states for t = ±oo, 



4 Of course the Schwinger-Dyson relation is also valid for the full propagator, including the border 
(13), but it is rather complicated to write it explicitly, and serves no great purpose here. 
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because there is no "vacuum" to the right, and hence the S-matrix formalism is not adequate 
in a world defined by such a Lagrangian. In order to avoid the latter, we can take a second 
border. In this case, one can easily see that the phase shift taken by the neutrino states keeps 
the S-matrix diagonal: the star is "transparent" for the neutrino propagation. Therefore, 
by applying eq. (pD|) , we can now conclude that the diagonality of the S-matrix justifies the 
null result for the weak energy density []. 

The generalization to the (3 + l)-dimensional problem is difficult even if the border is 
simplified to a flat one because of the presence of a refraction index, of a consequent modi- 
fication of the probability density of the waves at the border, including non-penetrating or 
non-outgoing waves. This (3 + l)-dimensional problem will be analytically and numerically 
studied in a forthcoming work [jToH . 

However, two main conclusions of this paper will remain valid in 3 + 1 dimensions: 
i) the natural connection between the neutrino sea and the border, ii) the proper definition 
of the ise infrared regulator for the propagator, or, equivalently, the correct definition of the 
zero energy level of the Dirac sea. On the other hand, the vanishing of the energy density 
w(z) will not remain valid in 3 + 1 dimensions. 

Finally, this work is a confirmation of the conclusion of ref. ||: after correcting the 
mistake of forgetting the pole in the analytic continuation, and adding the condensate, 
which is proved here to be directly connected to the border, the multibody exchange of 
massless neutrinos results in an infrared-well -behaved contribution for the neutron star 
self-energy 
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5 Despite the major analytic complexity introduced by the second border, we have computed the 
Feynman propagator and verified, by following the first method presented in this work, that the 
weak energy-density result is exactly zero. 
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Details of integration in 1 + 1 dimensions 
To compute the Feynman neutrino propagator in the momentum space, we should apply 
the results for the Fourier-transformed eigenstates (|T^) in eq. (|TT|). After some tedious 
transformations, we obtain: 



'■i 



j f + OO 

— / dk (h + I 2 + I 3 + h) , (21) 

Z7T JO 



where 



j = l x 

\(qi-ko)(qi-k +b)+ie 



i u,o / , / l^ „,o i 



|(7°-7 1 )7 



Uq%-k ){qi-k +b)-ie 2 W 1)1 q -k +K 



w- 



qi+ko)[qi+ko— bj+ie 
-\{ql+k )(ql+k -b)-ie 2 W ^ I J I q -k +ie ' 

3 ~~ \{qi+k Q )[qi+k +b)+ie 

, I I f 7 ° - 7 1 ) 7° I 

\{ll+ko){ql+k +b)-ie 2 V ' ' / ' qo+fc -Je ' 



(22) 



-■^(i?z-fco)(i?z-feo-fc)+i£ 
-|(9j-fco)(?|-feo-6)-ie 2 V / T I ) I q +ko-ie ' 

By taking into account the form of the different expressions for Jj in eqs. (p2[), we have: 

/ dk (h + I 3 ) = / dk h , 

JO JCi 

/•OO /■ 

/ d/c (J 2 + h)= dk I 2 , (23) 

JO JCi 



where the integration in the r.h.s should be done in the complex plane by following the 
contour C\ (see fig. H). 

Identifying the relevant poles for each case and applying Cauchy's theorem, we obtain 
the 1 + 1 propagator given above by eq. (|13|). 
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Contour Ci: ImCM 

(a) q >0 




(b) q o < 



Imfko) 
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FIG. 2. Integration contour C\. (a) qo > case; (b) go < case. 
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